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The elastic scattering cross section of ^^Mo(a,Q:)^^Mo has been measured at energies of i?c.m. ~ 
13, 16, and 19MeV in a wide angular range. The real and imaginary parts of the optical potential 
for the system ^^Mo - q have been derived at energies around and below the Coulomb barrier. The 
result fits into the systematic behavior of a-nucleus folding potentials. The astrophysically relevant 
^^Ru(7,a)^^Mo reaction rates at Tg — 2.0 and Tg — 3.0 could be determined to an accuracy of about 
16 % and are compared to previously published theoretical rates. 
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I. INTRODUCTION 



The nucleosynthesis of nuclei above the iron peak {A « 
60) proceeds mainly by neutron capture in the so-called 
s- and r-processes. In principle, neutron-deficient nuclei 
in this mass region (the so-called p nuclei, see Ref. (1) for 
a complete list) can be produced from more neutron-rich 
seed nuclei either by the removal of neutrons or by the 
addition of protons. However, due to the Coulomb bar- 
rier proton capture is strongly suppressed. There is gen- 
eral agreement that heavy neutron-deficient nuclei with 
masses above A « 100 have been synthesized by pho- 
todisintegration of previously produced nuclides at suffi- 
ciently high temperatures of (2 - 3) x lO^ K (Tg = 2 - 3, 
with Tg being the temperature in billion degrees). This 
so-called 7-process or p-process is discussed in detail in 
|0-|l2| . Several astrophysical sites for the 7-process have 
been proposed, whereby the oxygen- and neon-rich lay- 
ers of type II supernovae seem to be good candidates 
|0,|o|, but also exploding carbon-oxygen white dwarfs 
have been suggested [0. However, no definite conclu- 
sions have been reached yet. 

Nucleosynthesis calculations for the 7-process require 
a huge number of reaction rates. Up to about 1000 nuclei 



and 10000 reaction rates have been included in previous 
reaction networks [||. Recently, the complete network of 
the first self-consistent study of the 7-process including 
all relevant nuclei up to Bi amounted to about 3000 nuclei 
and all their respective reactions ||l^,|l^. Unfortunately, 
almost none of these reaction rates have been measured 
and the astrophysical calculations have to rely completely 
on statistical model calculations (e.g. Refs. [p^-p^). Of 
special importance are (7,Q!)/(7,n) branchings which de- 
termine abundance ratios of certain nuclides which, in 
turn, can in some cases be compared to abundances 
found in meteoritic inclusions p"8|-p0[| . It has been stated 
that the uncertainties for (7, li) reaction rates are huge 



||16,^ 21 1 . The determination of the a-nucleus poten- 
tial at astrophysically relevant energies helps to reduce 
the uncertainties of the calculations significantly |^,^ . 
(7,n) reaction rates have been measured in a recent work 
using a quasi-thermal photon spectrum, and rough agree- 
ment between theoretical predictions and the measured 
rates was found (2|j2^ . 

The overall agreement between the calculated and the 
observed abundance patterns of the p nuclei is relatively 
good. However, the mass region 70 < A < 100 is 
generally underproduced in the nucleosynthesis calcula- 
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tions The production among others depends 

on the neutron-producing ^^Ne(a,n)^^Mg reaction rate 
(which rnay enhance the s-process seed nuclei for the 7- 
process [12] ) and on the photodisintegration rates in the 
7-process but it remains unclear whether the underpro- 
duction can be cured by a change in those rates P6[ . 
Other explanations for this discrepancy include addi- 
tional production mechanisms like neutrino-induced nu- 
cleosynthesis Q and additional production sites like the 
rapid proton capture (rp) process on accreting neutron 
stars (e.g. |28 2^). However, it is still an open question 
whether rp-material can be ejected into the interstellar 
medium in sufficient quantities from these X-ray bursters 

The motivation of the experimental determination of 
the a-nucleus potential for ^^Mo is twofold. The de- 
termination of the a-nucleus potential at energies be- 
low the Coulomb barrier is limited in general because (i) 
the experimental data show only small deviations from 
the Rutherford cross section, and (ii) the optical poten- 
tials have ambiguities. An experiment on ^^Mo allows 
to extend the systematic study of a-nucleus potentials 
p2| , p3 31 3^ to lower energies. The second motivation 
refers directly to the production of the p nuclei ^^Mo and 
^■*Mo. A possible reaction path leading to the produc- 
tion of ^^Mo and ^''Mo is shown in Fig. |l]. Photodisin- 
tegration reactions of the nucleus ^^Ru can lead to the 
production of {i) ^^Mo by ^^Ru{-/ ,af^Mo and (ii) ^^Mo 
by 96Ru(7,n)95Ru(7,n)94Ru(2 x f3+ f^Mo. If this reac- 
tion path were the only production mechanism for ^^Mo 
and ^"^Mo, the abundance ratio between ^^Mo and ^^Mo 
would be directly related to the ratio of (7,0?) and (7,n) 
reaction rates of ^^Ru. In this case the ratio ^''Mo/^^Mo 
could be a thermometer for the 7-process because of the 
temperature dependence of the (7,n)/(7,a) branching ra- 
tio. However, for a quantitative analysis contributions of 
the rp-process to ^^Mo and ^"^Mo and the weak s-process 
contribution to ^"^Mo have to be known. 

The choice of the measured energies at about 13, 
16, and 19MeV has the following reasons. The astro- 
physically relevant energy window for (7,0;) reactions at 
Tg « 2 - 3 is of the order of « 8 - 10 MeV correspond- 
ing to 6 — 8 MeV for the reverse reaction ^^Mo(q;,7)®^Ru. 
Scattering experiments at these low energies are possible; 
however, a reliable determination of optical potentials is 
impossible because of the dominating Coulomb interac- 
tion. The height of the Coulomb barrier is about 15 MeV. 
We decided to measure at several energies above and be- 
low the Coulomb barrier to extract the optical potential 
and its energy dependence at energies as close as possible 
to the astrophysically relevant energy range. 

In the following paper we first present our experimental 
setup (Sect. ^). The experimental results are analyzed 
using systematic folding potentials, and discrete and con- 
tinuous ambiguities are discussed in detail (Sect. HI). 
The optical potential at astrophysically relevant energies 
is determined by extrapolation using the systematic be- 
havior of a-nucleus potentials ]3l]j3^, and the (7, a) re- 



action rates are calculatedfSect. [V). Finally, some con- 
clusions are given (Sect. |V|). A preliminary analysis of 
this experiment was presented already in [B3[. 



II. EXPERIMENTAL SETUP AND PROCEDURE 



The experiment was performed at the cyclotron lab- 
oratory at ATOMKI, Debrecen. We used the 78.8 cm 
diameter scattering chamber which is described in de- 
tail in Ref. |34|. Here we discuss only those properties 
which are important for the present experiment. A simi- 
lar setup has been used in our previous ^^*Sm(Q;,Q;)^ 
experiment [E2|. 



A. Targets 



The ^^Mo targets were produced by evaporation of 
97.33% enriched ^^MoOa on a thin (20/ig/cm2) carbon 
backing directly before the beamtime at the target labo- 
ratory of ATOMKI. The target was mounted on the tar- 
get holder in the center of the scattering chamber. The 
surface of the evaporated ^^MoOs turned out to be not 
very flat leading to relatively broad low-energy tails in 
the spectra of the elastically scattered a particles (see 
Fig. H). The target which was used during the whole 
experiment had a thickness of about 200/xg/cm^. The 
target stability was monitored during the whole experi- 
ment to avoid systematic uncertainties from changes in 
the target foil. 



B. Scattering Chamber 



A remote-controlled target ladder was placed in the 
center of the scattering chamber. Additionally, two aper- 
tures were mounted on the target holder to check the 
beam position and the size of the beam spot directly at 
the position of the target. The two apertures had a width 
and height of 2x6 mm^ and 6x6 mm^, respectively. The 
apertures were placed at the target position instead of 
the ^^MoOs target before and after each variation of the 
beam energy and beam current. The beam was opti- 
mized until no current could be measured on the larger 
aperture, and the current on the smaller aperture was 
minimized (typically less than 1 nA compared to about 
300 nA beam current) . The width of the beam spot was 
smaller than 2 mm during the whole experiment, which 
is very important for the precise determination of the 
scattering angle. Note that the relatively poor determi- 
nation of the height of the beam spot does not distu rb th e 
claimed precision of the scattering angle (see Sect. |IID|). 
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C. Detectors and Data Acquisition 

For the measurement of the angular distribution we 
used four siheon surface-barrier detectors with an active 
area A = 50 mm^ and thicknesses between D = 300 /im 
and D = 1500 /zm. The detectors were mounted on 
an upper and a lower turntable, which can be moved 
independently. On each turntable two detectors were 
mounted at an angular distance of 10°. Directly in front 
of the detectors apertures were placed with the dimen- 
sions 1.25 mm x 5.0 mm (lower detectors) and 1.0 mm 
X 6.0mm (upper detectors). Together with the distance 
from the center of the scattering chamber d = 195.6 mm 
(lower detectors) and d = 196.7 mm (upper detectors) 
this results in solid angles from Afi — 1.63 x 10^"* to 
AJ7 = 1.55 X 10"'''. The ratios of the solid angles of the 
different detectors were determined by overlap measure- 
ments with an accuracy better than 1 %. 

Additionally, two detectors were mounted at the wall 
of the scattering chamber at a distance of d = 351.3 mm 
and at a fixed angle of i? = 15° (left and right side rel- 
ative to the beam direction). These detectors were used 
as monitor detectors to normalize the measured angular 
distribution and to determine the precise position of the 
beam spot on the target. The solid angle of both monitor 
detectors is Ail = 8.10 x 10"*^. 

The signals from all detectors were processed using 
charge-sensitive preamplifiers (PA), which were mounted 
directly at the scattering chamber. The output signal 
was further amplified by a main amplifier (MA) and fed 
into ADCs. The data were collected using the commer- 
cially available system WinTMCA which provides an au- 
tomatic deadtime control which was found to be reliable 
in a previous experiment For the coincidence mea- 
surements (Sect. |lD) additionally the bipolar signals of 
the MAs were fed into Timing Single Channel Analyz- 
ers (TSCA), and the unipolar outputs of the MAs were 
gated using Linear Gate Stretchers (LGS). 

The energy resolution of the detectors was tested be- 
fore the experiment using a mixed a source and values 
better than 20keV were measured. During the experi- 
ment the achieved energy resolution is determined mainly 
by the energy spread of the primary beam and by the 
thickness and flatness of the target. Depending on the 
measured angle, the achieved energy resolution was be- 
tween 0.5 % and 2 % corresponding to AE « 200 keV at 
Ea ~ 20 MeV. Two typical spectra are shown in Fig. |^. 
The relevant peaks from elastic ^^Mo-a scattering are 
well separated from inelastic and from background peaks. 



D. Angular Calibration 

The angular calibration of the setup is of crucial impor- 
tance for the precision of a scattering experiment at en- 
ergies close to the Coulomb barrier because the Ruther- 
ford cross section depends sensitively on the angle with 



sin * {'d/2). A small uncertainty of 0.1° in the determi- 
nation of ■& leads to a cross section uncertainty of 2.0% 
(1.0%; 0.6%) at an angle d = 20° (40°; 60°). The fol- 
lowing methods were applied to measure the precise scat- 
tering angle {). 

The position of the beam on the target was continu- 
ously controlled by two monitor detectors. The precise 
position of the beam spot was derived from the ratio of 
the count rates in both monitor detectors. Typical cor- 
rections were smaller than one millimeter, leading to cor- 
rections in d of the order of 0.1°. However, because of the 
minor beam quality at the 16 MeV measurement, larger 
corrections had to be applied to this angular distribution 
leading to larger uncertainties in the determination of the 
optical potential. 

The position of the four detectors was calibrated using 
the steep kinematics of ^H(a,a)^H scattering at forward 
angles (10° < d\sh < 15°) iQ. The results of our previ- 
ous experiment could be confirmed within the uncertain- 
ties 11 . 

Finally, we measured a kinematic coincidence between 
elastically scattered a particles and the corresponding 
^^C recoil nuclei using a carbon backing without molyb- 
denum as target. One detector was placed at 'diah.a — 70° 
(left side relative to the beam axis), and the signals from 
elastically scattered a particles on ^^C were selected by 
a TSCA. This TSCA output was used as gate for the sig- 
nals from another detector which was moved around the 
corresponding ^^C recoil angle t^iab, recoil = 45.5° (right 
side). The maximum recoil count rate was found almost 
exactly at the expected angle (see Fig. ||). 

In summary, the overall uncertainty of the angles t9 in 
this experiment is about 0.1° for the measurements at 13 
and 19 MeV and about 0.2° for the 16 MeV measurement. 



E. Experimental Procedure and Data Analysis 

Three angular distributions have been measured at en- 
ergies of Ea = 19.50, 16.42, and 13.83 MeV. The beam 
current was between 80 nA and 320 nA. The experiment 
covers the full angular range from forward angles of 
"i? — 20° to backward angles of-d — 170° in steps of about 
1° at all three energies. The statistical uncertainties of 
each data point vary from < 0.1% at forward angles to 
about 1 % — 2 % at backward angles. 

The count rates iV(t9) in the four detectors have been 
normalized to the number of counts in the monitor de- 
tectors Nuoni^ = 15°): 
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with An being the solid angles of the detectors. These 
measured cross sections have been transferred to the 
center-of-mass system. The cross section at the moni- 
tor position i?Mon = 15° is given by the Rutherford cross 
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section. The relative measurement eliminates the typi- 
cal uncertainties of absolute measurements which come 
mainly from changes in the target and from the beam 
current integration. Nevertheless, the beam current was 
measured by standard current integration in the Fara- 
day cup, and the absolute value of the cross section was 
consistent with the measured relative cross sections. 

The three angular distributions are shown in Fig. ^. 
The lines are the result of optical model calculations (see 
Sect. [II). The measured cross sections cover five or- 
ders of magnitude between the highest (forward angles at 
E = 13MeV) and the smallest cross section (backward 
angles a.t E = 19MeV) with almost the same accuracy. 
Further details of the experimental set-up and the data 
analysis can be found in Ref. 36 1. 



III. OPTICAL MODEL ANALYSIS 

The theoretical analysis of the scattering data was per- 
formed in the framework of the optical model (OM). The 
complex optical potential is given by 



U{r) = Vc(r) + V{r) + iW{r) 



(3.1) 



where Vc{r) is the Coulomb potential, and V{r) and 
W{r) are the real and the imaginary part of the nuclear 
potential, respectively. 



A. The Folding Potential 



assuming identical proton and neutron distributions. For 
N Ri Z nuclei up to ^°Zr {Z ^ 40, N ^ 50) this assump- 
tion worked well |3^. However, to take the possibility 
into account that the proton and neutron distributions 
are not identical in the nucleus ^^Mo {Z ^ 42, N ~ 50) 
a scaling parameter w for the width of the potential has 
been introduced, which remains very close to unity. 

For a comparison of different potentials we use the inte- 
gral parameters volume integral per interacting nucleon 
pair Jfi and the root-mean-square (rms) radius ri.,„s,R, 
which are given by: 



Jr 



1 



ApAt 



V{r) d^r 



JV{r) r2 d^r 
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1/2 



(3.4) 
(3.5) 



for the real part of the potential V{r), and correspond- 
ing equations hold for W{r). Ap and At are the nu- 
cleon numbers of projectile and target. Note that in 
the discussion of volume integrals J usually the nega- 
tive sign is neglected; also in this paper all J values are 
actually negative. The values for the folding potential 
Vi (with A = w = 1) are Jfl = 267.88 MeVfm^ and 
rrms.R = 4.989 fm. 

The Coulomb potential is taken in the usual form of 
a homogeneously charged sphere where the Coulomb ra- 
dius i?c is chosen identically with the rms radius of the 
folding potential Vf. Rc = n-ms.R = 4.989 fm. 



The real part of the optical potential was calculated 
by a double-folding procedure: 



B. The Imaginary Potential 



Vtir) 



Different parametrizations of the imaginary part of 
pp{rp) pri^T) Vcs{E, p — pp + pT,s — |r + rp -tlKf |^pft?iKg,](ifsoitential were chosen. Volume and surface 

Woods-Saxon (WS) potentials are defined by the foUow- 
(3.2) ing equations 



where pp, pT are the densities of projectile and tar- 
get, respectively, and VcS is the effective nucleon- 
nucleon interaction taken in the well-established DDM3Y 
parametrization . Details about the folding proce- 

dure can be found in Refs. |^,^. The folding integral 
in Eq. ( ^^ ) was calculated using the code DFOLD ]4^ . 

The resulting real part of the optical potential V{r) 
is derived from the folding potential V{{r) by two minor 
modifications: 



V{r) = A • Vf{r/w) 



(3.3) 



Firstly, the strength of the folding potential is adjusted 
by the usual strength parameter A with A w 1.1 — 1.3. 
This leads to volume integrals of the real potential [see 
Eq. (|^)] of about Jr « 320 - 350MeVfm3 in the ana- 
lyzed energy range (22|P,||. Secondly, the densities of 
the a particle and the '''•^Mo nucleus were derived from the 
experimentally known charge density distributions pi] . 



with 



Wy{r) = Wo-f{x) (3.6) 
Ws{r)^Wo-^f{x) (3.7) 



/(x) = (l + e^)-i with x = (r-i?)/a (3.8) 



The depth Wq, the radius parameter R, and the diffuse- 
ness a have been adjusted to the experimental elastic 
scattering data. 

Fourier-Bessel (FB) potentials are given by 

n 

WpBir) = X! {knr/RFB)/{kTTr/RFB) (3.9) 

k=l 

with the cutoff radius -Rfb- Again, the Fourier-Bessel 
coefficients ak are adjusted to the experimental data. 
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C. Results and Continuous and Discrete Ambiguities 

1. The 19MeV Data 

The elastic scattering cross sections were calculated 
from optical potentials with the computer code AO p2[ . 
The code allows a variation of the potential parameters 
and determines the best-fit values from a test. 

In the first analysis the potential parameters A and w 
of the real part were kept close to the expected values 
from the systematic study | |2^ , |3l| , p^ : A « 1.1 — 1.3 and 
w w 1.0. Several parametrizations of the imaginary po- 
tential were tested. It was found that different imaginary 
potentials reproduced the experimental data with simi- 
lar quality. Five different fits are shown in Fig. || and 
the potential parameters are given in Tables |, || and p| . 
It turns out that the real potential is practically identi- 
cal in all these fits, but the shape of the imaginary part 
shows strong variations. The five imaginary potentials 
are shown in Fig. ^. 

Since we want to determine the optical potential at 
astrophysically relevant energies we have to extrapolate 
from the present measurements. Because of the oscillat- 
ing behavior of the Fourier-Bessel potentials we decided 
to use the combination of a volume and surface Woods- 
Saxon potential as basis for the extrapolation. The 
obtained with this potential is practically identical to the 
obtained from the Fourier-Bessel potentials. The cal- 
culations with a pure volume Woods-Saxon or a pure 
surface Woods-Saxon show significantly worse x^ values. 

In a second analysis the strength parameter A and the 
width w of the real potential were varied in a wider range. 
It was found that it is not possible to get a good fit to the 
data when the w parameter deviates from 1.0 by more 
than a few per cent. However, a variation of A leads 
to the known so-called "family problem". It is possible 
to obtain comparable fits to the experimental data with 
various A parameters. This phenomenon was discussed 
in detail for a similar experiment in [^ . In Fig. ^ we 
present the x^ values which were obtained from the fol- 
lowing procedure: (i) the parameter A was varied from 
about 0.5 to 3.5; {ii) the width parameter w and the 
imaginary part of the potential (consisting of a combina- 
tion of volume and surface Woods-Saxon potentials) were 
adjusted to the experimental data for each value of the 
strength parameter A. One can clearly see the families 2, 
3, 4, and 5 as minima in x^, corresponding to A values of 
0.81 — 1.52. Note that the minima are more pronounced 
than in the previous ^^'^Sm(a,Q;)^^'*Sm experiment p^ ] 
because the ratio E/Vc between the energy E and the 
Coulomb barrier Vc is much higher in this ^^Mo experi- 
ment. 

It is not possible to extract the correct family from 
these experimental data only. But together with the 
systematic behavior of the volume integrals found in 
|2|,|l]j32| we can decide that "family no. 4" (A = 1.256) 
should be used for the description of the experimental 



scattering data and for the extrapolation to astrophys- 
ically relevant energies (see Sect. iV). As mentioned 
above, family no. 4 with Jn « 340MeVfm'^ corresponds 
to the values of about Jr k, 320 — 350MeVfm^ which 
are expected from the systematics of a-nucleus poten- 
tials and also from other systems p3|-^5|. Neither fam- 
ily no. 3 with Jn « 280MeVfm'^ nor family no. 5 with 
J/j w 400MeVfm'^ fit into the systematics. Numerical 
problems in the fitting routine showed up for very shallow 
and very deep real potentials, and a clear determination 
of families 1 and 6—10 was not possible. 

One further interesting fact has to be mentioned. The 
real potentials corresponding to the families 1 — 10 are 
shown in Fig. |[ The potentials from families 2 — 5 which 
are well-defined as minima in x^ (see Fig. ^ have the 
same depth V{r) — — 2.66MeV at the radius r = 8.52 fm. 
However, not all potentials which have this depth do 
describe the data equally well; additionally, one has to 
find a minimum in x^ in Fig- 0- This behavior of a so- 
called "one-point potential" has been observed in sev- 
eral experiments, and the relevant radius has also been 
called "sensitivity radius" (see e.g. |46)); however, to our 
knowledge the additional restriction of a significant min- 
imum in x^ has only been observed in the analysis of the 
^'^''Sm(Q!,a)^'^^Sm data so far |^ which has been per- 
formed at a similar energy. 



2. The 13 and 16 MeV Data 



The procedure described in the previous Sect. HI CI 



was repeated for the lower energies of 13 and 16 MeV, 
and similar results were obtained. The imaginary part of 
the potential for all energies is taken as a combination of 
volume and surface Woods-Saxon potentials. The poten- 
tial parameters are listed in Table | and the calculations 
have been compared to the experimental data already in 
Fig. |4[ However, the 16 MeV data do not fit very well into 
the systematic behavior shown in Fig. ^. The potential 
extracted from the 16 MeV data has larger uncertainties 
than at the other energies because of experimental prob- 
lems (see Sect. II D). 



3. 30 Me V Data from Literature 

Two angular distributions at energies of about 30 MeV 
are available in literature 48|. Unfortunately, both 



angular distributions show systematic deviations between 
each other, and both distributions have not been mea- 
sured in the full angular range, but in the ranges 15° < 
i9c.m. < 75° @ and 15° < z9c.m. < 95° [Q. If one ad- 
justs the potential parameters to these discrepant angular 
distributions, discrepant optical potentials are obtained. 
The potential parameters are labeled in Table IV. How- 
ever, the potentials extracted from the data of |47|] fit into 
the systematics, whereas the data of Ea] do not fit. In 
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both cases the hmited angular range restricts the sensi- 
tivity of the potential parameters significantly. In Fig. ^ 
only the volume integrals derived from [^Tf are shown. 



4- Backward Angle Excitation Function 

The excitation function for ^^Mo(a,Q;)^^Mo has been 
measured by Eisen et al. at ??iab = 170° from 7 to 
16 MeV. We have calculated this excitation function from 
our best-fit potentials at 13, 16, and 19 MeV, and we find 
excellent agreement between the experimental and the 
calculated excitation function. The measured and the 
calculated excitation function at i^iab = 170° are shown 
in Fig. |l^. The calculation with the 13 MeV potential un- 
derestimates the deviation from the Rutherford cross sec- 
tion at higher energies. Such a behavior can be expected 
because of the smaller volume integral of the imaginary 
part in the 13 MeV data corresponding to weaker absorp- 
tion. However, the measured scattering cross sections at 
one special backward angle do not contain enough infor- 
mation to fix the optical potential and its energy depen- 
dence. 



5. (a,n)-induced reactions 



A set of experimental data corresponding to the re- 
action 92Mo(a,n)95Ru is available in |§ (Sect.IV B, 
Fig. 6a). We have calculated the cross section from our 
model in the measured energy range, and found very good 
conformity between our calculations and the existing ex- 
perimental data. However, the available data from the 
different experiments show discrepancies, which make it 
difficult to fix the energy dependence of the optical po- 
tential. Also the existing reaction data do not cover the 
astrophysically interesting energy range (between 7 and 
9 MeV, see Tab. VH), which would be helpful in order to 
confirm our predictions (see Sect. IV). 



jRiEc.m.) = Jr,o X exp[-(£;c 



Eo)yA^] (3.10) 



D. Discussion 

Volume integrals for various a-nucleus potentials in a 
broad range of masses and energies are shown in Fig. ^ 
for the real (|^A) and the imaginary part (^B and ^p) of 
the optical potential. The systematic behavior of volume 
integrals is also confirmed for lighter target nuclei [ ^j5l| 
and in various other systems which have been analyzed 
recently p3|-p5[|. For the extrapolation of the optical 
potential to astrophysically relevant energies (Sect. IV) 
parametrizations of the real and imaginary volume inte- 
grals are needed. 

As can be seen from Fig. there is only a weak 
energy dependence of the real volume integral Jj^ at en- 
ergies below the Coulomb barrier. As well as in Ref. [ p2| , 
a Gaussian parametrization is adjusted to the new data 



with JK,o = 337MeVfm3, £;o = 21.55 MeV and 
A = 147.01 MeV, leading to a curve (full line) which 
is somewhat flatter than the one proposed in Ref. p^ ] 
(dotted line). The uncertainties for extrapolations to 
lower energies are of the order of less than 5 % corre- 
sponding to about 10 — 20MeVfm'^. Furthermore, the 
shape of the real p otential is given by the folding pro- 
cedure (Sect. HI A). This means that the real part of 
the a-nucleus optical potential can be determined at en- 
ergies below the Coulomb barrier with relatively small 
uncertainties because {i) continuous ambiguities can be 
avoided using the folding potential and (m) discrete am- 
biguities can be resolved from the systematic behavior of 
a-nucleus potentials. 

The situation for the imaginary part of the potential 
is much worse. The volume integral J/ of the imagi- 
nary part depends strongly on the energy because many 
reaction channels open at energies around the Coulomb 
barrier. Different parametrizations have been proposed 
p0| , p2 53 1 . As an example we present the so-called 
Brown-Rho (BR) parametrization [p2[ 







Jo 



(Bo 



for 

for Er 



>Eo 



(3.11) 



with the excitation energy Eq of the first excited state. 
The saturation parameter Jq and the rise parameter A 
are adjusted to the experimentally derived values. An- 
other Fermi-like parametrization of the imaginary volume 
integral, first introduced in Ref. [E0|, reads 



JliEc.m.) 



Jo 



l + c^p[iE* - Ec.n,.)/a*] 



(3.12) 



with a similar saturation value Jq and the parameters E* 
and a* . The latter shape was also used for an attempt to 
derive a global a-potential ||5^ , with E* and a* depend- 
ing on Eq. However, the line derived with the parameters 
given in shows clear deviations from the new ^^Mo 
data. Therefore, we have adjusted this Fermi-like func- 
tion to the experimental data. Both parametrizations 
utilizing our fit parameters are shown in Fig. ^ (Brown- 
Rho) and (Fermi-like) for our new ^^Mo data. The 
parameters are listed in Table In the following, we 
will always use the parameters given in that table for the 
two descriptions unless specified otherwise. 

In general, the shapes of the BR and the Fermi 
parametrizations are quite similar: there is a saturation 
value Jo and a parameter that describes the steep rise 
of J/: A for BR and a* for the Fermi shape. How- 
ever, there are also important differences because the 
BR parametrization leads to a somewhat flatter rise of 
Jj than that of the Fermi function, and consequently, 
the extrapolation to lower energies is lower for the Fermi 
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parametrization than for BR. Con seque nces of these 
smaU differences will be given in Sect. IV C . Nevertheless, 
it should be noted that the BR parametrization only con- 
tains two free parameters, because Eq is fixed, whereas 
the Fermi shape has three free parameters, in principle. 

The shown ambiguities do not allow to determine the 
shape of the imaginary part. These ambiguities reduce 
the reliability of extrapolations to lower energies. A more 
stringent determination of the shape of the imaginary 
part requires extremely precise scattering data over a 
wide range of energies. A scattering experiment at about 
50MeV might help to reduce these uncertainties and to 
find the best parametrization of the imaginary volume 
integrals. 



IV. EXTRAPOLATION TO ASTROPHYSICALLY 
RELEVANT ENERGIES 

A. The Astrophysically Relevant Energy for (7,0) 
reactions 

The astrophysical decay rate is given by 



Jo 



(4.1) 



with the speed of light c, the cross section (T(^ Q,)(i?) of 
the (7,0;) reaction, and the photon density n^{E, T) of a 
thermal photon bath at temperature T 



n^{E,T) 



E^ 



he) exp(£;/fcr)-i 



(4.2) 



The integrand of Eq. (4.1) can be analyzed under the 
assumption that the astrophysical S-factor of the reverse 
(a, 7) reaction is constant: S(a.^^{E ) — const. Then the 
maximum of the integrand in Eq. (4.1) is found at the 
energy 



i;o(7,a) =-Bthr + ^G^' ( ^ 



kT 



2/3 



(4.3) 



with Eg — 2ii{'KZpZTe'^ /TiY and the threshold energy 
-Ethr for the (7,0;) reaction. The most effective energy 
i?o(7, oi) for (7,0;) reactions is given by the energy of the 
well-known Gamow window Eo{a, 7) for the inverse (a, 7) 
reaction plus the separation energy i^thr of the a particle. 
Note that the energy Eo{j,a) is the energy of the pho- 
ton, whereas the energy i?o(a,7) is the centcr-of-mass 
energy in the system ^^Mo - a. The astrophysically rele- 
vant energies for the system ^^Mo - a are listed in Table 

vii. 



In all astrophysical applications reaction rates are in- 
put only for reactions with positive Q value and the in- 
verse rate is then computed by applying detailed balance 
(see e.g. |15 ). That way, numerical stability of the re- 



of forward and reverse rates can be attained for a given 
channel. The rates in the two directions depend linearly 
on each other and thus the change of, say, the a potential 
equally influences both, in our case the a capture as well 
as the photodisintegration with a emission. This relation 
is valid provided that stellar rates are used in both direc- 
tions, accounting for thermal excitation of the respective 
targets. Because of that fact, in the following sections we 
make use of rate ratios so that the conclusions apply to 
the forward and inverse rates as well. 



B. Extrapolation of the Optical Potential 



As stated in Sect. HID, the extrapolation of the real 
potential can be performed reliably leading to Jn w 
325MeVfm'^ at astrophysically relevant energies with an 
uncertainty of about 5%. The corresponding strength 
parameter is A « 1.2. The width parameter was fixed at 
w = 1.0. 

The extrapolation of the imaginary part was per- 
formed as follows. In a first step the volume integral 
Ji was determined from the BR and Fermi parametriza- 
tions leading to J/ = 23.9MeVfm3 (BR) and J/ = 
15.4MeVfm^ (Fermi) at -Ec.m. = 5.81 MeV (correspond- 
ing to Tg — 2.0). The average of these values is Jj — 
19.6 ± 4.2MeVfm'^ which was used for the following cal- 
culations. 

The shape of the potential was taken as sum of volume 
and surface Woods-Saxon potentials where the radius pa- 
rameter R and the diffuseness a were estimated from the 
experimental data. The contribution of the volume term 
to Ji is assumed to be 30%, and the surface term con- 
tributes to 70%. This ratio is determined from the ex- 
perimental scattering data a,t E = 13, 19, and 30 MeV. 
The effect of a variation of the relative contributions of 
volume and surface term to J/ will be discussed below. 
These and other variations of the potentials allow an es- 
timation of the uncertainties of the calculated reaction 
rates. 



C. The ^''Ru(7,q)^^Mo reaction rate 

The variation of the reacti on ra tes when using various 
potentials is shown in Tables VIII and IX. In Table VIII 



action network is guaranteed and the proper equilibria 



the ratios of rates obtained with the different potentials 
in respect to a standard rate (taken from Refs. (l^Jl^ and 
using an a potential from Rcf. |^J]) are shown. As can 
be seen, the rate calculated with the global potential of 
Ref. [53) is lower by about two orders of magnitude than 
the standard rate. However, it was already stated above 
that this potential does not describe the ^^Mo data at 
higher energies. A simple equivalent square well poten- 
tial also yields a factor of 8 lower cross sections but 
neither does it describe the data nor is it considered to 
be reliable for this application |26| . When using the two 



7 



potentials with the extrapolated parameters for Tg = 2.0 
and Tg = 3.0 from Table 0, a reduction of the rate of 
about 40 - 50 % is found. 

Case AB explores the dependence on the geometry of 
the potential. A change in the geometry parameters of 
only 0.1 — 0.7% leads to a variation in the ratio of 7 — 10 % 
in the rate ratios which underlines the importance of ad- 
ditional scattering experiments to determine the shape 
of the imaginary optical potential. However, it should 
be mentioned that case AB is not fully consistent within 
our approach because it has a slightly different volume 
integral J/ and rms radius due to the unchanged depths 
of the volume and surface terms, but the differences are 
only of the same order of magnitude as those in the ge- 
ometry parameters. 

The sensitivity of the rates to variations in the extrap- 
olated volume integral J/ and the relative contributions 
of volume and surface term are studied in Tabic [X . Here 



the varied rates are compared to the rate obtained in case 
A of Table VHI . The ratios are given in the temperature 
range 0.5 < Tg < 10.0 in order to show the temperature 
dependence of those effects although strictly speaking the 
potential was derived assuming Tg = 2.0. 

The contribution of the surface term to J/ was varied 
within a reasonable range of 70 ± 20 %. This resulted in a 
variation of the rate of about ±6 %. Another uncertainty 
is introduced by the fact that we assumed the extrapo- 
lated Ji to be the mean between the value obtained by 
the BR and Fermi parametrizations. Using the higher 
BR value of J/ = 23.9 MeVfm'^ increases the rate by 
8% while using the lower value of Jj — 15.4 MeVfm'^ a 
suppression by about 10% is obtained. Thus, the error 
introduced by the different shapes of the parametriza- 
tions used for the extrapolation of Jj to low energies is 
dominating but still within satisfactory accuracy. 

Closing this section we conclude that the recommended 
rates are case A for Tg = 2.0 and case B for Tg = 3.0 from 



Table VIII with an error of 16%, mainly introduced by 



the ambiguities of the extrapolation of the imaginary part 
down to the relevant energies. The recommended rate is 
roughly a factor of two lower than the standard rate given 
in previous tabulations ||lq,0| . 



V. SUMMARY 



We measured the elastic scattering cross section of 
^^Mo(q;,q;)^^Mo in a wide angular range at energies of 
Ec.m. ~ 13, 16, and 19MeV. Additionally, data from lit- 
erature have been analyzed [^Tj-Q. The real and imag- 
inary parts of the optical potential for the system ^^Mo 
- a have been extracted from the data and extrapolated 
down to the astrophysically relevant energies around and 
below the Coulomb barrier. The result fits well into the 
known systematic behavior of a-nucleus folding poten- 
tials. The extrapolation of the imaginary part is not 
unique but our study shows that the use of two different 



energy dependencies introduces an error in the obtained 
rate of not more than 15 %. 

The derived stellar rates (for ^^Mo(q!,7)^^Ru as well as 
^^Ru(7,q;)^^Mo) are 50 — 60 % of the rates given in Refs. 
| p5| , p7t at stellar temperatures Tg — 2.0 — 3.0. Assuming 
the ™Ru(7,n)^'^Ru rate to remain unchanged, this would 
lead to a corresponding decrease in the abundance ratio 
^32]y[o/y94]y[o wlth rcspcct to an abundance ratio calcu- 
lated with the previous rate in the 7-process (as e.g. in 
|]l3| , p^ ). It is interesting to note that many network cal- 
culations for the 7-process [0-0 show an overproduction 
of ^^Mo relative to ^''Mo which may be reduced by the 
results of this work. However, as mentioned in the intro- 
duction, a complete analysis has not only to follow the 
7-process consistently but also to account for the possi- 
ble rp- and s-process contributions. This is beyond the 
scope of this paper. 
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FIG. 1. Nucleosynthesis path for ^^Mo and ''"Mo in the 
astrophysical 7-process. Stable nuclei are gray shaded. The 
nuclei ^2 Mo and ^"^Ru are marked; in this paper we determine 
an improved reaction rate for the ^®Ru(7,a)^^Mo reaction. 
Note that the (7,n) reactions stop at the magic neutron num- 
ber A'^ = 50. The circles at the right end of the diagram mark 
several other stable nuclei. 
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FIG. 2. Typical spectra zX. E ^ 19MeV and i9„ = 25° 
(upper diagram) and ~ 13 MeV and = 145° (lower 
diagram). The peak from elastic ^^Mo-a scattering is well 
separated from ^^C-a, elastic scattering (upper) and 

from inelastic ^'^Mo - a scattering (lower diagram). 
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FIG. 3. Relative yield of ^'^C recoil nuclei in coincidence 
with elastically scattered a particles. The shaded area shows 
the angle and the uncertainty which is expected from the cal- 
ibration using the steep kinematics of ^H(a,a)^H. The dotted 
line is a Gaussian fit to the experimental data points to guide 
the eye. 
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FIG. 4. Experimental cross sections of ''^Mo(a,Qf)''^Mo 
at iJc.m. ~ 13, 16, and 19 MeV normalized to the Ruther- 
ford cross section. The lines are the result of optical model 
calculations (see table ^. 
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5. Calculated cross sections of ^■^Mo(a,Q:)^'^Mo at 
Ec.m. ~ 19MeV using five different parametrizations of the 
imaginary part of the potential. These five fits look very sim- 
ilar; however, fits 1, 4, and 5 have a significantly smaller 
(see Tables | and |l|). 
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FIG. 9. Volume integrals of the real (A, upper) and 
imaginary part (B and C, center and lower diagram) of 
the optical potential derived from ^^Mo(Q,a)^^Mo scattering. 
For comparison volume integrals derived from ^°Zr(Q,Q)^°Zr, 
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from a emitters j32| were added. The lines in the up- 
per diagram show Gaussian parametrizations of the new 
^'^lAo{a,af'^Mo data (full line) and from Ref. Q (dotted 
line). The lines in the center diagram show the results of BR 
parametrizations of the imaginary part. In the lower di- 
agram the lines are the result of a Fermi-like parametrization 
from Grama and Goriely (GG) , using the parameters pro- 
posed by the authors, and the ones derived in this work (full 
line). Details see text. 
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TABLE I. Potential parameters of the imaginary part of 
the optical potential (combination of volume and surface WS 
parametrizations) , derived from the angular distribution of 
^^Mo(q,q)^^Mo at E =13, 16 and 19MeV, and its integral 
potential parameters J and rrms of their real and imaginary 
parts. 

fit E Wy (MeV) Ry (fm) ay (fm) Wp (MeV) Ro (fm) ao (fm) 



1 


19 


-1.584 


1.7667 


0.2659 


334.25 


1.2605 


0.2073 




16 


-9.558 


1.668 


0.248 


308.20 


1.348 


0.099 




13 


-5.128 


1.656 


0.002 


467.06 


1.369 


0.071 



fit E 


A 


w 


Jr 


"^rms.R 


Ji 






(MeV) 






(MeVfm^) 


(fm) 


(MeVfm^) 


(fm) 




1 19 


1.257 


1.003 


337.2 


4.991 


86.2 


5.806 


2.15 


16 


1.346 


0.9974 


357.8 


4.976 


85.9 


5.992 


4.84 


13 


1.352 


0.9758 


336.7 


4.869 


67.3 


6.043 


1.26 



TABLE II. Potential parameters of the imaginary 
part of the optical potential (volume and surface WS 
parametrization) , derived from the angular distribution of 
^^Mo(a,a)^^Mo at E =19 MeV, and its integral potential pa- 
rameters J and rrms of their real and imaginary parts. 

fit Wy (McV) Ry (fm) ay (fm) Wp (MeV) Rp ( fm) ap (fm) 

2 -10.806 1.7116 0.3601 

3 121.76 1.4947 0.4206 



fit A W Jr rrms.R J I Trms,! ~)d I F 

No. (MeVfm^) (fm) (McVfm") (fm ) 

~2 1.237 1.010 341.1 5.037 57^9 6l31 3.67 

3 1.188 1.021 338.9 5.095 80.6 6.957 4.28 



TABLE III. Potential parameters of the imaginary part of 
the optical potential (FB para- metrization) , derived from the 
angular distribution of ^^Mo(a,Q:)®^Mo at E =19 MeV, and 
its integral potential parameters J and rrms of their real and 
imaginary parts. 

fit -Rfb (fm) ffll 02 03 04 05 06 07 

~^ 12li 154.94 -311.98 839.41 -325.97 880.41 - 
5 12.0 167.94 -329.11 1103.43 -592.10 1666.61 -346.79 1118.06 



fit 


A 


w 


Jr 


Trms.R 


Ji 






No. 






(MeVfm=^) 


(fm) 


(MeVfm^^) 


(fm) 




4 


1.272 


0.998 


338.8 


4.979 


68.2 


4.524 


2.23 


5 


1.287 


0.991 


336.0 


4.947 


53.9 


4.319 


2.14 
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TABLE IV. Potential parameters of the imaginary part 
of the optical potential (combination of volume and surface 
WS parametrizations) , derived from the angular distribution 
of ^'^Mo{a,af^Mo from Refs. Q and Q, and its integral 
potential parameters J and rrms of their real and imaginary 
parts. 

Ref. Wv (MeV) Rv (fm) ay (fm) Wo (MeV) Ro (fm) ao (fm) 



[47 -4.91 1.78 0.40 183.13 1.13 0.36 

[48 -3.51 1.72 0.95 337.06 1.26 0.23 



Ref X W Jr 7-rms,R Jl rmis,! 

No. (MeVfm^) (fm) (MeVfm^) (fm) 



[47 1.19 1.022 334.76 5.098 90.95 5.704 
[48 1.15 1.014 315.49 5.056 107.91 6.014 



TABLE V. Extrapolated values of the potential paramters 
of the imaginary part of the optical potential at the astro- 
physically relevant energies Eo — 5.8 MeV (Tg = 2.0) and 
Eo — 7.6 MeV {Tg = 3.0), and its integral potential parame- 
ters J and Trnis of their real and imaginary parts. 



E Wv (MeV) Rv (fm) ay (fm) 
7.6 -1.466 1.717 0.268 
5.8 -1.091 1.720 0.270 



Wo (MeV) R o (fm) ao (fm) 
36.86 1.295 0.419 
27.44 1.297 0.420 



E 


A 


w 


Jr 




Ji 




(MeV) 






(MeVfm^) 


(fm) 


(MeVfm^) 


(fm) 


7.6 


1.219 


1.000 


327.1 


4.991 


26.2 


6.085 


5.8 


1.209 


1.000 


324.3 


4.991 


19.6 


6.095 



TABLE VL Parameters of the BR and Fermi parametriza- 
tions of the imaginary volume integral Ji for ^■^Mo(Q,a)^^Mo. 

parametrization saturation value rise parameter other parameters 
BR Jo = 99.8 MeVfm^ A = 7.68 MeV Eo = 1.51 MeV 

Fermi Jo = 91.0 MeV fm^ a* = 2.78 MeV E* = 10.24 MeV 



TABLE VIL Most effective energies Eo for the 
Mo(a,7)''^Ru and the ^'^Ru(7,a)^^Mo reactions. 



To 


Eo(a,f) 


Eo{l,a) 


2.0 


5.81 


7.51 


2.5 


6.75 


8.44 


3.0 


7.62 


9.31 
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TABLE VIII. Ratio ^ = r^/ro of the astrophysical reac- 
tion rates obtained with different imaginary potentials to 
a standard rate ro ]l5| , p^ . 

Tg Gg^ ESW*^ A (r9 = 2) B (Tg = 3) AB ° 
2.0 0.014 0.121 0.453 0.497 0.503 

3.0 0.012 0.140 0.546 0.579 0.585 

''Potential from Ref. H] 

'^Equivalent square well potential, e.g. as in Ref. 
'^Imaginary part from the extrapolated values at i5 = 5.8 MeV 
in Tab. 0. 

''imaginary part from the extrapolated values at _E = 7.6 MeV 
in Tab. 0. 

''Potential depths are from A, the geometry parameters from 
B. 



TABLE IX. Variation of the imaginary part of the poten- 
tial derived for Tg = 2.0. Rate ratios are shown in respect 
to the rate obtained with the parameters for E = 5.8 MeV 
{J I = 19.6 MeVfm^, 70% surface contribution) from Tab. 

Variation of Ji Surface Contribution 



T 


J I = 


J I = 


J I -- 


= 19.6 


MeVfm'^' 


10^ K 15.4 MeVfm^ 


23.9 MeVfm^ 


90% 


80% 


60% 


50% 


0.5 


0.942 


1.077 


0.952 


0.976 


1.019 


1.043 


1.0 


0.949 


1.070 


0.949 


0.975 


1.032 


1.057 


1.5 


0.913 


1.078 


0.937 


0.968 


1.029 


1.057 


2.0 


0.902 


1.077 


0.937 


0.969 


1.029 


1.056 


2.5 


0.918 


1.062 


0.945 


0.974 


1.026 


1.049 


3.0 


0.945 


1.050 


0.958 


0.981 


1.020 


1.040 


3.5 


0.968 


1.027 


0.966 


0.984 


1.016 


1.030 


4.0 


0.991 


1.013 


0.977 


0.988 


1.010 


1.021 


4.5 


1.007 


1.001 


0.983 


0.992 


1.006 


1.014 


5.0 


1.019 


0.995 


0.989 


0.995 


1.005 


1.009 


6.0 


1.028 


0.991 


1.000 


1.000 


1.000 


1.009 


7.0 


1.031 


0.985 


0.999 


0.999 


1.000 


1.000 


8.0 


1.026 


0.987 


1.000 


1.000 


1.000 


0.996 


9.0 


1.020 


0.990 


1.004 


1.002 


0.998 


0.996 


10.0 


1.011 


0.994 


1.005 


1.002 


0.998 


0.994 
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